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Abstract. 

A comprehensive study of relativistic and resonance effects in electron impact 
excitation of (e+Fe XVII) is carried out using the Breit-Pauli R-matrix (BPRM) 
method in the relativistic close coupling (CC) approximation. Two sets of 
eigenfunction expansions are employed; first, up to the n = 3 complex corresponding 
37 fine-structure levels (37CC) from 21 LS terms; second, up to the n= 4 corresponding 
to 89 fine-structure levels (89CC) from 49 LS terms. In contrast to previous works, 
the 37CC and the 89CC collision strengths exhibit considerable differences. Denser 
and broader resonances due to n= 4 are present in the 89CC results both above 
and below the 37 thresholds, thus significantly affecting the collision strengths for 
the primary X-ray and EUV transitions within the first 37 n= 3 levels. Extensive 
study of other effects on the collision strengths is also reported: (i) electric and 
magnetic multipole transitions E1,E2,E3 and M1,M2, (ii) J-partial wave convergence 
of dipole and non-dipole transitions, (iii) high energy behaviour compared to other 
approximations. Theortical results are benchmarked against experiments to resolve 
longstanding discrepancies — collision strengths for the three prominent X-ray lines 
3C, 3D and 3E at 15.014, 15.265, and 15.456 A are in good agreement with two 
independent measurements on Electron-Beam- Ion- Traps (EBIT). Finally, line ratios 
from a collisional-radiative model using the new collisional rates are compared with 
observations from stellar coronae and EBITs to illustrate potential applications in 
laboratory and astrophysical plasmas. 
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1. Introduction 

The ground configuration of Fe XVII has a stable closed L-shell structure (neon core) 
rendering Fe XVII the dominant Fe ion species in many laboratory and cosmic plasmas. 
Owing to its importance Fe XVII has also become a benchmark ion, and one of the 
most extensively studied for its spectral diagnostic potential in plasma- and astro- 
physics. Both the atomic structure of Fe XVII and electron impact excitation (EIE) of 
(e + Fe XVII) are extremely complex. Although EIE of Fe XVII has been investigated 
experimentally and theoretically for a long time Q, there are considerable uncertainties 
in line intensity comparisons between measurements and theories, even with the most 
elaborate methods. 

Unlike lower ionization stages of Fe, Fe XVII is a medium-to-highly charged 
many-electron ion with strong correlation and relativistic effects in both the target 
dynamics and the electron collision processes @, particularly for non-dipole forbidden 
and intercombination transitions and near-threshold cross sections. Both intermdediate 
coupling LSJ, and jj-coupling, are needed in atomic structure calculations for proper 
target level designations. Precise atomic structure calculation is a crucial part of EIE 
calculations. However, despite a number of experimental and theoretical works, accurate 
atomic structure information on Fe XVII is not available in literature. The first 27 levels 
of Fe XVII in the configurations 2p 5 31 (l=s,p,d) have been determined experimentally. 
But consistent energies for inner-shell excitation configurations 2s2p 6 31 (l=s,p,d) (from 
level 28 to 37) are yet to be obtained. The same is true for the n= 4 configurations 
(2p 5 41 (l=s,p,d,f); 2s2p 6 41 (l=s,p,d,f)). Systematic theoretical structure calculations 
have been used to determine energies not available from experimental measurements, 
with an estimated accuracy of less than a few to ten per cent, and an uncertainty of 
~0.lA in wavelengths to the ground state. 

Transition probabilities and electron impact excitation collision strengths of 
Fe XVII are more inaccurate. We calculate transition probabilities using the code 
SUPERSTRUCTURE , which provides the structure input for the EIE calculations of 
Fe XVII || , and the multi-configuration Dirac-Fock (MCDF) method using the GRASP 
code with extensive configuration-expansion Q. All El, Ml, E2 transitions among the 
89 levels are calculated using SUPERSTRUCTURE and GRASP codes. We find some 
large M2 and E3 A-values from the comprehensive GRASP data, and for the first time, 
point out their importance in spectral formation of Fe XVII . 

Similar to the study of the Fe XVII structure, there is a long history of EIE 
calculations for Fe XVII . Basically, there are three types of EIE calculations in literature, 
(i) Distorted- wave (DW) and relativistic distorted- wave (RDW), which couple initial 
and final channels for a given transition and thus only background cross sections are 
calculated, (ii) Isolated resonance approximation plus (R)DW: including a limited 



4 



number of resonances in the isolated resonance approximation (IRA). This approach 
improves (i) somewhat; however we show that the final results are not very accurate 
due to the inevitably limited number of channels in the IRA, and breakdown of IRA 
due to overlapping of dense and broad resonances in Fe XVII |J. In some cases, the 
results from this approach may be misleading or confusing owing to unknown missing 
channels, (iii) Close coupling (CC) approach: non- relativist ic R-matrix calculations in 
LS coupling, and relativistic BPRM calculations in intermediate coupling. While this is 
the most advanced method, and in principle capable of including all important atomic 
effects, several issues require extreme care, as demonstrated in this work. Interestingly, 
there is no complete R-matrix calculations in literature for Fe XVII either. One previous 
coupled-channel calculation by Mohan etal. |§ employed 15 LS terms, followed by a pair 
coupling transformation to obtain fine-structure collision strengths between levels in the 
first 27 levels. However, Mohan etal. obtained only the background collision strengths, 
similar to the various DW calculations, because their calculations were above the highest 
threshold and therefore no resonances were included. Nontheless their data have been 
used to compare with experiments and observations, and the conclusion that resonance 
enhancement is not important has been drawn. Recently, limited BPRM calculations 
for Fe XVII were reported [[/J. Similar to the work in |J, only 27CC was employed 
in the calculations. Also, no detailed resonance structures in collision strengths were 
presented in this work either and thus it prevents the in-depth investigations. In order 
to address these and other issues thoroughly, it is therefore necessary to carry out a 
full set of BPRM calculations for (e + Fe XVII) . There are some R-matrix calculations 
for the neon isoelectronic sequence for other ions. Ne-like selenium has been calculated 
by fully relativistic Dirac R-matrix approach M and BPRM approach |J using a 27CC 
expansion in jj-coupling and in LSJ coupling, respectively. 

In the present work, relativistic and resonance effects are considered with two sets 
of calculations: 21 terms or 37 fine-structure levels (37CC) up to the n= 3 complex, and 
49 terms or 89 fine-structure levels (89CC) including up to the n— 4 levels. The 37CC 
BPRM calculation includes more than 7300 channels, with the largest Hamiltonian 
matrix of dimension 3659; the 89CC BPRM calculation involves more than 20,000 
channels overall, with and the largest the Hamiltonian matrix of dimension 10086 
(this is possibly the largest BPRM calculation yet carried out). There are pronounced 
differences between the two sets of calculations. The results show that the 89CC 
calculation with levels up to n= 4 is necessary to obtain accurate collision strengths 
even for the transitions up to the n= 3, since the n— 4 resonances appear both above 
and, surprisingly, below the n= 3 levels and affect the effective collision strengths 
significantly. The two sets of calculations also show that the backgrounds for some 
important transitions are also affected. The implications of this study on Fe XVII for 
Ne-like ions and other highly charged ions are pointed out. 
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As mentioned, in addition to theoretical studies there have been several experimental 
measurements for Fe XVII . For example, experimental line intensity ratios have recently 
been measured on Electron-Beam-Ion- Traps (EBIT) at Lawrence Livermore National 
Laboratory (LLNL, Brown etal. 1998), and at the National Institute for Standards 
and Technology (NIST, Laming etal. 2000), for some prominent and strong lines of 
Fe XVII that have been observed in the X-ray spectra of solar corona and other 
stellar coronae, active galactic nuclei, X-ray binaries, supernovae, and recently in solar- 
type star Capella from X-ray satellites Chandra and XMM-Newton [TO]. These 3 
lines correspond to (see Fig. 1) X-ray lines from M-shell 3d-2p with wavelength at 
~15l: (1) 3C A 15.0131: ls 2 2s 2 2p 5 [l/2]3d 3/2 *P° (level 27) — ► ls 2 2s 2 2p 6 % (ground 
state); (2) 3D A 15.2651: ls 2 2s 2 2p 5 [3/2]3d 5/2 3 D° (level 23) — ► ls 2 2s 2 2p 6 %; (3) 3E 
A 15.4561: ls 2 2s 2 2p 5 [3/2]3d 5/2 3 P? (level 17) — > ls 2 2s 2 2p 6 %. The line intensities 
display subtle effects since the 3C is a dipole-allowed transition, but the 3D and 3E lines 
are intercombination transitions that behave as forbidden transitions (spin forbidden) at 
low-Z neon-like ions, and as allowed (in jj coupling) El transitions for high-Z high ions. 
Comparison of experimental and theoretical data shows disagreement of up to 50% for 
R1=3C/3D, and a factor of two for R2=3E/3C JTl], 12|. This means that some atomic 
mechanisms, and their effect on line formation, have not been considered in previous 
theoretical works. 

The outline of the paper is as follows. In Sees. 2 and 3 the basic theoretical and 
computational methods and techniques are briefly described. A schematic illustration 
of many important Fe XVII lines are shown in Fig. 1, and discussed in Sec. 4. In Sees. 5 
and 6 we present the results for cross sections and line intensities, and demonstrate 
that: (i) dense resonance structure appeared in all transitions over the entire energy 
range below the highest threshold in the 89CC BPRM collision strengths; in particular, 
resonance enhancement generally dominates forbidden and intercombination transitions 
(but has not heretofore been studied), (ii) the theoretical line intensity ratios for 
intercombination transitions, over the strongest dipole-allowed transition, agree with two 
sets of recent EBIT measurements [11, [12[ to 10% or within experimental uncertainties. 
The conclusions are summarised in Sec. 7. 



2. Theory 

The Fe XVII target wavefunctions are computed using SUPERSTRUCTURE (Eissner 
etal. 1974 |3|]), which employs a scaled Thomas-Fermi-Dirac-Amaldi potential (Eissner 
and Nussbaumer 1969 [D3]) to compute the set of one-electron orbitals. The scaling 
parameters are optimised with a list of target configurations and the Breit-Pauli 
Hamiltonian [|TJj] . In the Breit-Pauli R- matrix (BPRM) approximation JT5| the following 
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Hamiltonian terms are retained 
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where H^ + i is the non-relativistic Hamiltonian together with the one-body mass 
correction term, the Darwin term and the spin-orbit term resulting from the reduction 
of the Dirac equation to the Pauli form. The mass-correction and Darwin terms do 
not break the LS symmetry, and they can therefore be retained with a great effect in 
computationally cheaper LS calculations. Spin-orbit interaction does, however, split 
the LS terms into fine-structure levels labelled by Jn, where J is the total angular 
momentum and 7r the parity. 

In the coupled channel or close coupling (CC) approximation the wave function 
expansion, ^(E), for a total spin and angular symmetry SLn or Jn, of the (N+l)- 
electron system is represented in terms of the target ion states as: 



where x% is the target ion wave function in a specific state SiL^i or level J(K^ and 9i is the 
wave function for the (N+l) i/l electron in a channel labeled as SiL^Jj)^ kfi^SLn) [Jit]; 
kf is the incident kinetic energy. In the second sum the <&_,-' s are correlation wave 
functions of the (N+l)-electron system that (a) compensate for the orthogonality 
conditions between the continuum and the bound orbitals, and (b) represent additional 
short-range correlations that are often of crucial importance in scattering and radiative 
CC calculations for each SLn. The $j's are also referred to as "bound channels", as 
opposed to the continuum or "free" channels in the first sum over the target states. In 
the relativistic BPRM calculations the set of SLir are recoupled to obtain (e + ion) 
states with total Jn, followed by diagonalization of the (N+l)-electron Hamiltonian. 
Details of the diagonalization and the R-matrix method are given in many previous 
works (e.g. Berrington etal. 1995). 

3. Computations 

3.1. Target eig en functions 

The configuration-expansion consists of 49 LS terms corresponding to 89 fine-structure 
levels with principal quantum number up to n= 4. Target energies are given in 
Table 1 and compared with observed ones. Based on both SUPERSTRUCTURE and 
GRASP calculations we have assigned both jj-coupling and intermediate coupling LSJ 
designations to the energy levels. The MCDF GRASP calculations are with the same 
configuration-expansion as used in the SUPERSTRUCTURE calculations. The MCDF 
method is a complete self-consistent-field (SCF) procedure, which means both the 
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orbitals and the expansion coefficients are variational. However, the Breit interaction 
is included as a perturbation in the Cl-type SUPERSTRUCTURE calculations with 
fixed orbitals. Selected oscillator strengths for the transitions among target states 
are compared with other works in Table 2. We also report on some higher multipole 
transition probabilities, compared to El, in Table 3. The extensive calculations on 
Fe XVII atomic structure will be reported separately. 

3.2. Electron impact excitation of Fe XVII 
Target configurations 

2s 2 2p 6 , 2s 2 2p 5 (3s,3p,3d), 2s 2 2p 5 (As, Ap, Ad, Af), 2s 1 2p 6 (3s, 3p, 3d), 2s 1 2p 6 (As,Ap,Ad,Af) 

corresponding to 89 fine-structure target levels up to n= A are included in the CC 
expansion. The largest symmetry is total J = 9/2 (even parity) with 400 free channels 
and 86 bound channels ($j); the dimension of the Hamiltonian is 10086 with 25 
continuum orbitals in the inner region. The R-matrix boundary is at Rq = 3.969 
a.u. The maximum electron impact energy is up to 400-500 Ry. BPRM collision 
strengths are calculated for all J < 51/2, including all the partial waves in the range 
I — — 31. This ensures convergence of the collision strengths without partial wave top- 
up in the low-energy region. In the high-energy region, we use relativistic distorted-wave 
(RDW), and/or Coulomb-Born-Bethe (CBe), top-up procedures for the different types of 
transitions: allowed, intercombination, forbidden, or other mixed or peculiar transitions. 
Special attention is paid to convergence with respect to resolution of resonances. We 
use a constant-energy of up to 20,000 energies to compute rate coefficients for practical 
applications. We also carry out a 37CC calculation with the n= 3 levels in order to 
check the accuracy of previous DW and other calculations, and to demonstrate the 
necessity of performing 89CC calculations. The other R-matrix parameters in the 37CC 
calculations were the same as the 89CC. 

4. Important Fe XVII lines 

To illustrate applications of the present Fe XVII results to astrophysical objects and laser 
transtiions in plasmas, we show selected X-ray transitions to the ground level, and X/UV- 
ray transitions between excited levels, in Figs. 1(a) and 1(b). We have computed 3653 
transition probabilities and collision strengths for 3916 transitions among the 89 levels 
for line intensity modeling in a collisional-radiative model (CRM). The 20 important 
lines are shown in Fig. 1(b). Some of them, dipole El transitions 3A (1-33) and 3B 
(1-31), E2S (1-37) transition from an upper level with an inner-shell 2s hole, the strong 
El line 3C (1-27), and two intercombination transitions 3D (1-23) and 3E (1-17) from 
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2p 5 3d levels to the ground state, show prominent emission and/or absorption features 
in collision ionized plasmas. These lines have been extensively used in the diagnostics 
of temperature, density and ionization balance, 

Fig. 1(b) also shows El allowed and intercombination lines 3G (1-3) and 3F (1-5) 
respectively, and a M2 (1-2) line, from low-lying upper levels of the 2p 5 3s configuration. 
The upper levels 3F and 3G have very fast radiative decay rates. This point is very 
important in the creation of population inversion of soft X-ray laser lines in Fe XVII 
. As opposed to the high-lying lines from 2p 5 3d levels and others, cascade effects may 
dominate the line formation of 3F, 3G, and M2. These three low-lying lines are powerful 
diagnostic tools used in photoionized plasmas. 

An E2L (1-7) line from 2p 5 3p has also been observed from EBIT experiments, and 
possibly from some cosmic plasmas. The X/UV line A 1153 A (3-4) between excited levels 
3 and 4 is a potential diagnostic of temperature, such as in solar flares, by comparing 
its intensity to the A 975 A line in Fe XVIII. 

There are two important points related to the X/UV-ray lines in Fig. 1(b) among 
excited levels from 2p 5 3p to 2p 5 3s: (1) fast radiative decay of levels 3 and 5, and 
(2) strong monopole collisional excitation rate to level 15 (2p 5 3p 1 S'o) responsible for 
population inversion leading to the strongest soft X-ray laser lines. In addition to the 
collisional excitation effects, cascade effects may also contribute to the level population 
of level 15. 

5. Results and discussions 



5.1. Fe XVII atomic structure 
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Table 1. The 89 fine-structure n = 2, 3 and 4 levels included in the BPRM calculation 
and their calculated and observed energies in Rydbergs for Fe XVII. 'obs' values 
are observed values from NIST website ( http:/ / www, nist. onj ) (level denoted by 'f 
is from Q. 'SS' and 'MCDF' are from SUPERSTRUCTURE and GRASP calculations, 
respectively. The index i is used for transition keys. NB:, both jj-coupling and LS- 
coupling notations are needed for the level designations of Fe XVII. 



i 


Level 




obs 


SS 


MCDF 


i 


Level 


obs 


SS 


MCDF 


1 


2s 2 2p 6 1 S (0,0)0 


0.0 


0.0 


0.0 


46 


2s 2 2p 5 4p 3 P 2 (3/2,3/2)2 




72.9701 


72.7671 


2 


2s 2 2p 5 3s 


3 P° (3/2 l/2l°2 


53.2965 


53.3542 


53.1684 


47 


2s 2 2p 5 4p 3 P (3/2,3/2)0 




73.2625 


73.0529 


3 


2s 2 2p 5 3s 


X P\ (3/2,l/2)°l 


53.43 


53.4982 


53.3100 


48 


2s 2 2p 5 4p 3 Di (1/2,1/2)1 




73.7627 


73.5563 


4 


2s 2 2p 5 3s 


3 P° fl/2 l/2)°0 


54.2268 


54.2786 


54.0957 


49 


2s 2 2p 5 4p 3 Pi (1/2,3/2)1 




73.8425 


73.6430 


5 


2s 2 2p 5 3s 


3 P° fl/2 1/2)°1 


54.3139 


54.3726 


54.1851 


50 


2s 2 2p 5 4p 1 P> 2 (1/2,3/2)2 




73.8597 


73.6595 


6 


2s 2 2p 5 3p 


3 Si (3/2,1/2)1 


55.5217 


55.5586 


55.3963 


51 


2s 2 2p 5 4d 3 P^ (3/2,3/2)°0 




74.0121 


73.8044 


7 


2s 2 2p 5 3p 


3 D 2 (3/2,1/2)2 


55.7787 


55.8321 


55.6606 


52 


2s 2 2p 5 4d 3 P° (3/2,3/2)°l 


73.95 


74.0462 


73.8369 


8 


2s 2 2p 5 3p 


3 D 3 (3/2,3/2)3 


55.8974 


55.9429 


55.7791 


53 


2s 2 2p 5 4d 3 F^ (3/2,5/2)°4 




74.0921 


73.8777 


9 


2s 2 2p 5 3p 


x Pi (3/2,3/2)1 


55.9804 


56.0297 


55.8654 


54 


2s 2 2p 5 4d 3 P° (3/2,5/2)°2 




74.1021 


73.8904 


10 


2s 2 2p 5 3p 


3 P2 (3/2,3/2)2 


56.1137 


56.1586 


56.9950 


55 


2s 2 2p 5 4p 1 S (l/2,l/2)°0 




74.0943 


73.9033 


11 


2s 2 2p 5 3p 


3 Po (3/2,3/2)0 


56.5155 


56.5799 


56.4050 


56 


2s 2 2p 5 4d 3 F° (3/2,3/2)°3 




74.1080 


73.8994 


12 


2s 2 2p 5 3p 


3 Di (1/2 1/2)1 


56.6672 


56.7209 


56.5495 


57 


2s 2 2p 5 4d X D° (3/2,3/2)°2 




74.1534 


73.9456 


13 


2s 2 2p 5 3p 


3 Pi (1/2,3/2)1 


56.9060 


56.9475 


56.7855 


58 


2s 2 2p 5 4d 3 D 3 (3/2,5/2)°3 




74.1864 


73.9736 


14 


2s 2 2p 5 3p 


X D 2 (1/2,3/2)2 


56.9336 


56.9798 


56.8135 


59 


2s 2 2p 5 4d 3 D° (3/2,5/2)°l 


74.30 


74.3841 


74.1666 


15 


2s 2 2p 5 3p 


1 So (1/2,1/2)0 


57.8894 


58.0832 


57.9308 


60 


2s 2 2p 5 4f 3 Di (3/2,5/2)1 




74.6730 


74.4521 


16 


2s 2 2p 5 3d 


3 P^ (3/2 3/2)°0 


58.8982 


58.9384 


58.7738 


61 


2s 2 2p 5 4f 1 G 4 (3/2,5/2)4 




74.6761 


74.4522 


17 


2s 2 2p 5 3d 


3 P° (3/2,3/2)°l 


58.981 


59.0169 


58.8454 


62 


2s 2 2p 5 4f 3 G 5 (3/2,7/2)5 




74.6799 


74.4535 


18 


2s 2 2p 5 3d 


3 PS (3/2 5/2)°2 


59.0976 


59.1644 


58.9826 


63 


2s 2 2p 5 4f 3 D 2 (3/2,7/2)2 




74.6855 


74.4623 


19 


2s 2 2p 5 3d 


3 F< (3/2 5/2)°4 


59.1041 


59.1820 


58.9901 


64 


2s 2 2p 5 4f 3 D 3 (3/2,7/2)3 




74.7109 


74.4871 


20 


2s 2 2p 5 3d 


3 Fg (3/2,3/2)°3 


59.1611 


59.2240 


59.0498 


65 


2s 2 2p 5 4f 3 F 2 (3/2,5/2)2 




74.7141 


74.4910 


21 


2s 2 2p 5 3d 


1 T>2 (3/2,3/2)°2 


59.2875 


59.3519 


59.1797 


66 


2s 2 2p 5 4f X F 3 (3/2,5/2)3 




74.7170 


74.4935 


22 


2s 2 2p 5 3d 


3 Dg (3/2,5/2)°3 


59.3665 


59.4492 


59.2598 


67 


2s 2 2p 5 4f 3 F 4 (3/2,7/2)4 




74.7246 


74.4992 


23 


2s 2 2p 5 3d 


3 D° (3/2,5/2)°l 


59.708 


59.7904 


59.6082 


68 


2s 2 2p 5 4d 3 F^ (l/2,3/2)°2 




75.0232 


74.8203 


24 


2s 2 2p 5 3d 


3 F2 (l/2,3/2)°2 


60.0876 


60.1439 


59.9749 


69 


2s 2 2p 5 4d 3 D° (l/2,5/2)°2 




75.0470 


74.8391 


25 


2s 2 2p 5 3d 


3 D 2 (l/2,5/2)°2 


60.1617 


60.2190 


60.0344 


70 


2s 2 2p 5 4d X F 3 (l/2,5/2)°3 




75.0704 


74.8618 


26 


2s 2 2p 5 3d 


1 F° (l/2,5/2)°3 


60.197 


60.2643 


60.0754 


71 


2s 2 2p 5 4d 1 PJ (l/2,3/2)°l 


75.17 


75.2381 


75.0263 


27 


2s 2 2p 5 3d 


1 P° (l/2,3/2)°l 


60.6903 


60.8342 


60.6279 


72 


2s 2 2p 5 4f 3 G 3 (1/2,5/2)3 




75.6139 


75.3975 


28 


2s2p 6 3s 3 


Si (1/2,1/2)1 




63.3223 


63.2125 


73 


2s 2 2p 5 4f 3 G 4 (1/2,7/2)4 




75.6249 


75.4057 


29 


2s2p 6 3s 1 


So (1/2,1/2)0 




63.7906 


63.6986 


74 


2s 2 2p 5 4f 3 F 3 (1/2,7/2)3 




75.6331 


75.4147 


30 


2s2p 6 3p 3 


P° (l/2,3/2)°0 




65.7252 


65.6346 


75 


2s 2 2p 5 4f 1 D 2 (1/2,5/2)2 




75.6316 


75.4155 


31 


2s2p 6 3p 3 


P° (1/2,1/2)°1 


65.601 


65.7608 


65.6676 


76 


2s2p 6 4s 3 Si (1/2,1/2)1 




81.6960 


81.5889 


32 


2s2p 6 3p 3 


P^ (l/2,3/2)°2 




65.9265 


65.8380 


77 


2s2p 6 4s 'So (1/2,1/2)0 




81.8528 


81.7414 


33 


2s2p 6 3p 1 


P? (l/2,3/2)°l 


65.923 


66.0704 


65.9782 


78 


2s2p 6 4p 3 P^ (l/2,l/2)°0 




82.6666 


82.5483 


34 


2s2p 6 3d 3 


Di (1/2,3/2)1 




69.0140 


68.9221 


79 


2s2p 6 4p 3 P° (1/2,1/2)°1 


82.52 


82.6786 


82.5594 


35 


2s2p 6 3d 3 


D 2 (1/2,3/2)2 




69.0332 


68.9323 


80 


2s2p 6 4p 3 P^ (l/2,3/2)°2 




82.7438 


82.6281 


36 


2s2p 6 3d 3 


D 3 (1/2,5/2)3 




69.0657 


68.9518 


81 


2s2p 6 4p X P° (l/2,3/2)°l 


82.67 


82.7908 


82.6724 


37 


2s2p 6 3d 1 


D 2 (1/2,5/2)2 


69.282f 


69.4386 


69.3247 


82 


2s2p 6 4d 3 Di (1/2,3/2)1 




83.9021 


83.7838 


38 


2s 2 2p 5 4s 


3 P^ (3/2,l/2)°2 




71.8688 


71.6517 


83 


2s2p 6 4d 3 D 2 (1/2,3/2)2 




83.9097 


83.7886 


39 


2s 2 2p 5 4s 


1 P° (3/2,l/2)°l 
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71.9139 


71.6983 


84 


2s2p 6 4d 3 D 3 (1/2,5/2)3 




83.9229 


83.7976 
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72.5803 
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84.0529 
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84.4858 


84.3522 


44 


2s 2 2p 5 4p 


3 D 3 (3/2,3/2)3 




72.8992 


72.6948 


89 


2s2p 6 4f X F 3 (l/2,7/2)°3 




84.4962 


84.3621 


45 


2s 2 2p 5 4p 


1 P 1 (3/2,3/2)1 




72.9287 


72.7243 
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Table 2. Comparisons of weighted oscillator strengths (gf) for selected dipole allowded 
transitions for Fe XVII in length- form (L) and velocity- form (V). 'SS' indicates present 
SUPERSTRUCTURE calculations; 'MCDF' results are obtained with GRASP; 'NIST' 
values are the recommended values from NIST website; 'CIV3' are from fl6| and 
'MBPT' are from 0; the quantity aeb stands for a x 10 b . 
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15 
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5.1.1. Energies The calculated energies in Table 1 from SUPERSTRUCTURE and 
GRASP are compared with observed values wherever available. Table 1 also provides 
the key to the level indices for transitions in tabulating transition probabilities, 
collision strengths, and Maxwellian averaged collision strengths. The comparison with 
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Table 3. Selected transition probabilities As of Fe XVII: the first and second entries 
for each El transition are the respective values of A in the length-form and velocity- 
form; the quantity aeb stands for a x 10 b . 
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experimental values shows agreement to within 1% or better. The target representation 
is thus generally better than in most EIE calculations. 

5.1.2. Transition probabilities We report here briefly our calcultions of transiton 
probabilities. Full results will be reported elsewhere due to the space limitation. Table 
2 presents weighted oscillator strengths for a number of El transitions, whose accuracy 
indicates that of the target eigenfunctions. Table 3 gives a few selected transition 
probabilities of dipole-allowed El transitions, forbidden magnetic dipole Ml and electric 
quadrupole E2 transitions, and some unexpectedly strong E3 and M2 transitions. Only 
E2 and E3 in length-from are tabulated in Tanle 3. In intermediate-coupling, or in jj- 
coupling, intercombination transitions are a special type of El transitions. The excellent 
agreement of the length-form (the first entry) and the velocity-form (the second entry) A- 
values for each El transitions selected in Table 2 is a further prove that we have obtained 
a good target for EIE calculations. In Table 3, all selected E3 and M2 transitions are of 
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order 10 4 s _1 or more. They may have some influence on the modeling of line emissions. 
In particular, the M2 line from level 2-1 has long been observed as a prominent line in 
astronomy and in laboratory photoionized plasmas, The population of level 2 is fed up 
by the cascade effects from 2p 5 3s, 2p 5 3p, and 2p 5 3d and other higher configurations. 
Accurate M2 transition probability is key to model this line. This line also has an 
important plasma diagnostics potential. 



5.2. Electron impact excitation of Fe XVII 

5.2.1. 37CC collision strengths In order to understand the dependence of resonances 
on the n-complexs, we first calculate a set of 37CC collision strengths. Also, there are 
many previous 37-level (R)DW calculations so it is more straightforward to compare 
previous results with the present 37CC calculations. 

Fig. 2 shows extensive resonances in several transitions. Fig. 2(a) is the collision 
strengths for the first resonant excitation Q( 1 S - 3 P 2 ), a magnetic quadrupole (M2) 
transition. The dominant role of resonances is clear, when compared with the 
background collision strengths and in particular compared to the earlier (R)DW 
calculations, which correspond only to the non- resonant background. Fig. 2(b) is for 
transition from the ground state to the 15th energy level 2p 5 3p 1 So- This transition is 
a monopole transition which is very important in soft X-ray lasing studies. The green 



dots are RDW calculations in [18|], and the green squares are our RDW calculations 



T9fl . Both RDW calculations agree very well with the background collision strengths of 



the BPRM results. 



5.2.2. Comparison with non-resonant approximations: We discuss in particular the 
RDW calculations as representative of a long list of previous calculations in literature 
using the distorted-wave (DW) method, the Coulomb-Born (CB) method, or their 
variants. Although the (R)DW values are in very good agreement with the background 
values in the present calculations, the dense resonance structures are not accounted for. 
However, the relativistic effects and correlation effects in the target are considered in the 
for in RDW calculations. The RDW method may be thought of as a 2CC calculation 
including initial and scattered channels, without the remaining channel coupling and 
resonance effects which, as we demonstrate in this work, are not only significant but 
may dominate certain transitions. Since DW results are still extensively used in spectral 
modeling, care should be exercised to take these limitations into account. 

The one previous non-relativistic R-matrix calculation also agrees well with the 
present background values in the high-energy region above the highest threshold in the 
27-level target. Collision strengths are reported at a few energies in Ref. in this 
region. It is clear that since all calculations are above the highest threshold in the CC 
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expansion, there are no resonance structures included in the previous R-matrix work. 

It should be noted that although for the 3 transitions 3C, 3D and 3E discussed in 
detail here, the RDW results are in excellent agreement with the background values 
of the 89CC BPRM calculations, there are very large discrepancies for some other 
transitions even in the background values due to broad diffuse resonance or background 
enhancement. The differences can be from a factor of two to up to an order of magnitude, 
depending on the strength of the transition, (see Figs. 12 for further discussions). 

5.2.3. Partial wave expansion: the 89CC collision strengths Before presenting the 
collision strengths from the 89CC calculations, it is instructive to point out the issues 
of the convergence of partial wave expansion addressed in detail in our study. We used 
the partial wave collision strength fij, and the partial sum of PW collision strengths 
Z)5y 2 ^J> as shown in Figs. 3,4 and 5. All partial cross sections are from the 89CC 
BPRM calculations unless otherwise indicated. 

In Figs. 3(a)-(d), Vtj and J2i/2^J f° r the transitions 3C and 3G are plotted as 
a function of 2J. For each transition, ttj and J2i/2 are calculated at four incident 
electron energies: 100 Ry (black curve), 200 Ry (green), 300 Ry (blue), and 400 Ry 
(Red). The filled circles in Figs. 3(b) and (d) are the corresponding total collision 
strengths that partial sum of PW collision strength should converge onto. It is 

noted that there is a peak and a minimum for all curves of flj for the 3C transition, and 
there are also some irregular features around 7 < 2J < 17. The higher partial waves 
(2 J > 53) not calculated in this work all fall in the tail part of the f2j curves. This is 
necessary to ensure accurate CBe top-up (see Sec. 3, also discussed later). It is found 
that for the 3C convergence is achieved; the largest difference is < 5% at = 400 Ry. 

Figs. 3(b) and (d) show the partial wave convergence for the 3F transition. Again, 
there is one minimum for most curves of flj; however, there are two peaks and two 
minima = 100 Ry. Furthermore, the 3F curves are much more irregular than the 
3C, as seen from Figs. 3. When we check the convergence from Figs. 3(b) and (d), we 
find that at low Ei convergence has been achieved. For higher Ei however, complete 
convergence is not obtained in the CC calculations. But the positions of the last peak 
for all the Vtj curves are well below the highest partial wave (2J=51) employed in the 
present 89CC calculations, so the CBe top-up approach should again be very accurate 
and we have convergent collision strengths for these types of transitions. 

The transitions shown in Figs. 3 are from the ground level to the n— 3 levels. 
In Figs. 4(a)-(d), we discuss further a particular transition, 2-44 (indices in Table 1,) 
between excited levels from n— 3 to n— 4. Similar features as the curves of Vtj above 
are also found in Fig. 4(a). Noticeable differences are: (i) the positions of the peaks are 
at much higher J values (close to the highest J used in CC calculations), in particular 
for high Ei energies; (ii) the values of fij at 2 J = 51 are still a big fraction of the peak 
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values. This feature implies that the CC results have certainly not converged, as can 
also be proved from the partial sum of J2{/2 m Fig- 4(b). In fact, at i?j=400 Ry, 
E1/2 is only ~50% of the converged value for the transition 2-44. The main reason 
for this difference, compared to Figs. 3, is that the transition energies for n— 3-4 are 
much smaller than for the n=2-3. Fortunately, the positions of all highest peaks for 
all ft j curves are located below (though close to) 2J=51, so the CBe top-up approach 
should be accurate. To reconfirm this, we did a RDW calculation for the transition 2-44 
including partial waves up to 80. The results are plotted in Figs. 4(c) and (d). We have 
a slightly different notation for the RDW results. In Fig. 4(c), the partial wave Qe is 
plotted as a function of the orbital angular momentum I of scattered free electron; while 
in Fig. 4(d), the partial sum of fli is plotted as a function I. However, according to the 
collision theory of partial wave expansion, we arrive at an identical conclusion. As seen 
from Fig. 4(c), convergent results have been obtained with partial waves up to £=80 for 
the transition 2-44. 

Partial wave analysis shown in Figs. 3-4 are for dipole allowed transitions. It is 
expected that the convergence for E2 transitions would also be very slow. This can 
be seen from the two E2 transitions, 2-51 and 6-60, shown in Figs. 5(a)-(d). Similar 
features to Figs. 4 are found, and the same conclusion can be drawn. In addition, we 
used the Burgess- Tully method |20| to study the convergence of the E2 cross sections 
in Figs. 6, where the reduced collision strengths Q r are plotted as a function of reduced 
electron energy E r = Ei/(Ei + c) (c=e is a constant, E r G [0,1]), for the E2 transitions 
2-51 and 6-60. The red filled circle corresponds to infinite incident energy; our RDW 
results are also shown for comparision as blue squares. From this figure we conclude that 
the BPRM calculations of E2 transitions 2-51 and 6-60 are well converged for Ej < 400 
Ry (E r <0.7). 



5.2.4- 89CC Collision strengths The 89CC calculations reveal the presence of 
resonances due to the n= 4 levels that appear not only above the n= 3 thresholds, 
but also below those energies. This in particular affects the comparison of theoretical 
results with the recent EBIT measurements of relative line ratios, as discussed in [|| . 

Below, we discuss several aspects of the large 89CC calculations. 

A. Comparison of 31 CC and 89CC results: 

The differences between the two sets of calculations are illustrated for the three 
transitions 3C, 3D, and 3E in Figs. 7-9, respectively, that correpsond to important X-ray 
lines observed in laboratory and astrophysical plasmas. The RDW calculations and non- 
relativistic R-matrix calculation || without resonances are also shown for comparison. 

The 3C is the strongest resonance transition from the ground state in the emission or 
absorption spectra of Fe XVII . Figs. 7(a) and 7(b) show the 37CC and the 89CC collision 
strength respectively. The green open square is the non-relativistic DW calculation [T[], 
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and the blue triangles are previous 27CC R- matrix calculations || [JJ. In Fig. 7(c), the 
37CC collision strengths (red) is overlapped with the 89CC results in order to show the 
difference between them. Three main differences are noticed: (i) Resonances in the range 
of 69 Ry<Ej <75 Ry show up in 89CC calculations due to the inclusion of targets states 
in configurations 2p 5 4£(£ =s,p,d,f). Practical implication of this point are discussed 
later, (ii) Denser resonances appear in the 89CC calculations for Ej <69 Ry. This is 
also important in comparing the relative line intensity measured in EBIT experiments, 
and in obtaining correct effective collision strengths under different electron velocity 
distributions (e.g. gaussian, maxwellian, numerical), (iii) The background collision 
strengths are lower in the 89CC calculations. This is due to the effect of channel 
coupling. Redistribution of electron fluxes at a larger number of thresholds, as opposed 
to a smaller number, reduces the background. Although there is only a few percent 
difference in the 3C background collision strengths, it is the strongest transition which 
is a factor of a few to orders of magnitude larger than other transitions. Therefore 
the redistribution of the 3C electron collision flux results in a significant, 10% or more, 
change for other transitions. 

The green and red dashed-lines shown in the plot are the numerically averaged 
(NA) and the gaussian averaged (GA with a 30 eV FWHM) collision strengths used to 
calculate line intensity ratios discussed in Sec. 6.1. The NA and GA values can also be 
used to assess resonance enhancements due to the 89CC calculations, from Figs. 7(a) 
and (b). 

Figs. 8 and 9 provide details of resonance enhancements for the intercombination 
transitions 3D and 3E, respectively. The RDW calculations |18[ and R-matrix 
calculations || 0] are also shown. Several features may be noted, (i) Resonance 
enhancements in 3D and 3E in the 89CC calculations are much more pronounced than 
for the 3C. Clearly, the n= 4 levels add considerably more resonances, both above and 
below the n= 3 complex. Again, similar to the case for 3C in Figs. 7, these effects 
are more readily discernible by comparing the NA and GA values from the respective 
37CC and 89CC calculations (relative to their respective backgrounds), (ii) While 
both the 3D and 3E are intercombination transitions, with denser resonances than the 
3C, the energy behavior of all three transitions is different. The 3D collision strength 
increases with Ej, with basically the same energy-dependence as the dipole allowed 3C, 
but the 3E decreases with increasing Ej over a very broad energy range, typical of a 
forbidden transition, (iii) The physically different energy behaviour is also related to 
the appropriate coupling scheme that should be used to describe these three transitions 
in each case. Whereas the dipole allowed transition 3C can be treated in either LS or 
jj-coupling scheme, the 3D is more appropriatly considered in a relativistic jj-coupling 
scheme. However, a non-relativistic LS-coupling scheme may describe the behaviour of 
the 3E (spin) forbidden transition. These effects demonstrate that the atomic structure 
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and collision dynamics in intermediate Z and medium-to-highly charged many-electron 
atomic species such as Fe XVII are rather complicated, since they lie in a transient 
region where relativisitc effects begin to manifest themselves. 

Finally, we note that the dominant role of the n— 4 resonances seen in the 
89CC results in Figs. 7-9, is due to resonant configurations 2s 2 2p 5 3£4£' that manifest 
themselves from ~47 Ry, considerably below the excitation thresholds of the 2p — > 
3d lines 3C, 3D, and 3E at ~60 Ry. Therefore there are relatively fewer n = 3 
resonances, and the n = 4 resonances greatly influence the near-threshold behavior 
of these cross sections. In the region 75 - 84.5 Ry the resonances are due to thresholds 
corresponding to two-electron-excitation configurations 2s2p 6 4£ with weakly coupled 
channels and resonances, and should be affected by the higher n > 4 resonances. 

B. Wavefuncion expansion coefficients: 

The ^-dependence of 3C, 3D, and 3E can also be quantitatively addressed from the 
eigen function expansion of the upper levels in the Fe XVII target. We have: upper level 
(27) of 3C: 0.7857|27) + 0.1753|23) + 0.0305|17); upper level (23) of 3D: 0.7479|23) + 
0.2010|27) + 0.0491|17); upper level (17) of 3E: 0.9150|17) + 0.0767|23) + 0.0030|27). 
The transitions probabilities for 3C, 3D and 3E are 2.47el3, 6.01el2, and 9.42el0, 
respectively. Level 17 is nearly described in pure LS-coupling from its eigen function 
expansion, while levels 23 and 27 begin to considerably depart from pure LS-coupling 
to intermediate-coupling scheme (or jj-coupling). It is the strong coupling of level 23 
with level 27 that makes the 3D line behave as a dipole-allowed transition, as evident 
from the mixing coefficients. In LS-coupling, the 3E is a spin-forbidden transition, so 
its collision strengths decreases with Ei. 

C. Other excited El transitions: 

In addition to the dipole allowed El transitions (including intercombination El 
transitions) discussed above, two other interesting El transitions 1-5 (3F) and 1-33 
(3A) are shown in Figs. 10. The upper level of the 3F (as the 3G in Fig. 3) is strongly 
mixed, therefore LS coupling designations are not appropriate, as discussed in Sec. 1. 
Instead, we need to use jj-coupling notation. Because of the strong coupling in the 
target, the two transitions show similar resonance structures. Both exhibit the allowed 
transition behavior, although the 3F is in fact an intercombination transition. The 
RDW results and previous R- matrix results are also shown for comparison (Fig. 10a), 
with good agreement for the background. 

The upper level of the transition 3A (Fig. 10b) is from a configuration with a 2s- hole, 
therefore the strong channel coupling effects with the 2s2p 5 4£ (£=s,p,d,f) configurations 
cause the full range of pronounced resonances up to the highest threshold of the n=4 
complex. 

D. Electric quadrupole (E2) and magnetic dipole (Ml) transitions: 

In Fig. 11, collision strengths for six forbidden transitions (pure E2, pure Ml, mixed 
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E2 and Ml) are displayed. In Figs. 11(a) and (b) two pure E2-type transitions from 
the ground level, 1-7 and 1-37, are shown. Resonant features in 1-37 are similar to 1-33 
since the upper level of 1-37 also has a 2s-hole (as discussed in subsection C above). In 
Fig. 11(a), one previous R-matrix value for 1-7 is smaller than the present background, 
and all RDW values, by more than a factor of two. In Figs. 11(c) and (d), two pure 
Ml-type transitions from the ground level, 1-6 and 1-9, are shown. In Fig. 11(d), one 
previous R-matrix value for 1-9 is larger than the present background and the RDW 
values by more than a factor of two. In Figs. 11(e) and (f), two mixed E2 and Ml-type 
transitions between excited levels 2-3 and 2-5, are shown. Although there is no 2s-hole in 
the upper levels 3 and 5, it is interesting that some prominent resonant features appear 
near the highest threshold of the n=4 complex. 

E. Electric octupole (E3) and magnetic quadruple (M2) type transitions: 

A particularly noteworthy finding of the present work is the potential significance 
of the E3 transitions in spectral modeling of Fe XVII based on their relatively large 
A- values (Table 3). In Fig. 12, four forbidden transitions (pure E3, pure M2, mixed 
E3 and M2, and mixed El, E3 and M2 types) are given. In Fig. 12(a), a pure M2- 
type transition from the ground level 1-2 is shown. This transition is also presented in 
Fig. 2(a) from the 37CC calculations. One previous non- resonant R-matrix value for 1-2 
is smaller than the present background and all the RDW values by more than a factor of 
two [F2l[] . The RDW values agree with the present background in the high-energy region. 
However, in the low-energy region, the RDW values differ from the present background 
by about a factor of two. In Fig. 12(b), a pure E3-type transition from the ground level 
1-20 is shown. All RDW results, and the previous non-resonant R-matrix values, are in 
good agreement with the present background. 

A mixed M2 and E3-type transition between excited levels 3-36 is shown in 
Fig. 12(c), and a mixed El, M2 and E3-type transition between excited levels 2-35 is 
shown in Fig. 12(d). As the cases discussed earlier, the upper levels of 3-36 in Fig. 12(c), 
and 2-35 in Fig. 12(d) have a 2s-hole leading to pronounced resonances in the energy 
range up to the n= 4 complex for both transitions. The RDW calculations are off by 
more than a factor of two for both transitions; while the previous DW values are differ 
by an order of magnitude, as shown in the figures. 

F. The Monopole transitions: 

Finally, we illustrate the monopole transitions of great interest in laser excitations 
in Fe XVII since they are collisionally, and not radiatively, excited. In Fig. 13, two 
monopole transitions are given. In Fig. 13(a), the monopole transition from the ground 
level 1-15 is shown. This transition is also presented in Fig. 2(b) from the 37CC 
calculations. In Fig. 13(b), the monopole transition from the ground level 1-11 is shown. 
The RDW values differ slightly but the previous non-resonant 27CC R-matrix values 
differ considerably with the present background collision strengths. 
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6. Rate coefficients and line intensity ratios 

Significant resonance enhancement of the collision strengths of forbidden and 
intercombination transitions has been demonstrated in this work. This directly enhances 
the rate coefficients which, in turn, affect the computation of line intensities. While the 
calculation of rate coefficients is a voluminous task, currently under way (to be reported 
another paper), we present selected results from this study and apply them to the 
analysis of previous experiments and observations. 



6.1. Comparisons with EBIT experiments 

As described earlier, two sets of EBIT experiments have measured the relative intensities 
of the 3C, 3D, and 3E transitions [II, 0. The monoenergetic beams used in these 



experiments sample the effective cross sections averaged over the beam width. Although 
the beam width is relatively large (>30 eV FWHM) compared to the resonance widths, 
unlike astrophysical plasmas with typically Maxwellian distribution over a wide range 
of electron velocities, the monoenergetic EBIT experiments accurately probe the major 
atomic processes in line formation: electron impact excitation and radiative decay, and 
are not affected by other processes. Line intensities measured in the EBIT experiments 



for 3C, 3D and 3E transitions can be modeled as 11, 22 



l kl 



rj(\ kl )B kl ((T kl ■ v)n e N Fe XV n (3) 



where the upper state k is 27, 23 and 17 respectively for the three transitions. Response 
function ^(A^i) of the spectrometer is taken to be the same for all the three lines because 
the wavelengths of the three lines — 15.0L4, 15.26 A, 15.46 A — are close together. 
The branching ratios for all the radiative decay routes are calculated to be 1.0, 1.0, 0.89 
for 3C, 3D and 3E, respectively. (cr fel t>) is the beam averaged rate coefficient; n e and 
Np e xvii are the electron density and the density of Fe XVII ions respectively. One 



then obtains 22 



Rl = 3C/3D = (<r 3C ■ v)/{<t 3D ■ v) (4) 
R2 = 3E/3C = 0.89{a 3E -v)/(a 3C -v) (5) 

For monoenergetic electron beams employed in EBIT experiments, we use two ways 
to calculate the rate coefficients using the collision strengths in Figs. 7-9 with complex 
resonance structures. The first is a direct numerical average (NA), and the second is 



a Gaussian average (GA) using a beam width of 30 eV JTT|. The averaged results are 
not very sensitive either way, owing to the fairly large beam widths, although there 
is a small variation in some regions. Both the NA and GA averaging results are also 
shown in panels (a)&(b) in Figs. 7-9 as red and green curves respectively. From the 
working assumptions outlined above, the ratios of line intensities in EBIT experiments 
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are in fact the ratio of their averaged-collision strengths convolved over the electron 
beam width, as in Eqs. (4) and (5). We note particularly that the observed up and 
down 'oscillation' in the line ratio Rl seen in the EBIT experiments [23| is likely not 
from the experimental error bars, but from the averaged-collision strengths with varying 
resonance enhancements in different energy regions, and/or different types of transitions. 
Both, the collision strengths for the 3C, 3D, or 3E, and the line ratios Rl and R2 as a 
function of electron energy, have additional complexities because of fine-structure and 
coupling effects. Therefore, constant or simple ratios of collision strengths at one energy, 
extrapolated to the entire energy range, are not accurate. 



Table 4. Comparison of present line ratios for R1=3C/3D and R2=3E/3C with EBIT 
measurements 



£i=0.85 keV 0.9 keV 1.15 keV 





EBIT 


2.77±0.19 a 


2.94±0.18 6 


(3.15±0.17,2.93±0.16) a 


R1=3C/3D 


Theory c :NA 


2.80 


3.16 


3.05 t 




Theory c :GA 


2.95 


3.27 


3. 10t 




Other Theory 


3.78 d ; 


4.28 e ; 


3.99 / 




EBIT 




0.10±0.01 6 




R2=3E/3C 


Theory :NA 


0.11 


0.085 


0.07 t 




Theory c :GA 


0.11 


0.083 


0.07 t 




Other Theory 


0.04 d ; 


0.05 e ; 


0.05 f 



a EBIT experiments at LLNL [|T 

d 



with NA and GA; 



1 



/ 



6 EBIT experiments at NIST ||12|| ; c present theory 
H; ^present values with extrapolation of resonance 



enhancement from ab initio collision strengths from E < 1.02 keV (see text). 



Table 4 shows the line intensity ratios calculated from Eqs. (4) and (5) and the 
procedures discussed above. The agreement with two independent EBIT experiments is 
excellent, to 10% or within the experimental error bars. In the transitions 3C, 3D, and 
3E we find strong resonances in the 89CC calculations for all energies up to 75.0 Ry (1020 
eV), corresponding to level 75 (Table 1), indicating that channel couplings in this region 
are accounted for. However, even higher thresholds (e.g. n— 5) are likely to contribute 
to resonances above this threshold. In order to compare with EBIT experimental data 
at these higher energies, we extrapolate the resonance enhancement to the region above 
75.0 Ry. The predicted values are also given in Table 4. It should be mentioned that 
although the highest threshold in the 89CC calculations is ~84.5 Ry, more practically 
we define the highest strong 'resonance' threshold to be at 75 Ry because the levels 
76-89 are from configurations 2s2p 6 4£ (£ = s,p,d,f) and the channel-couplings with 
transitions 3C, 3D, and 3E is therefore very weak, as expected from Figs. 7-9 in the 
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energy range E t > 75 Ry |2"4]| . The physical implication of these weak couplings is that 
the flux transfer between channels arising from two-electron-excitation processes is very 
weak in the lines 3C, 3D, and 3E. 

6.2. Comparison with observations from stellar coronae 

Next, we carry out the line intensity calculation using a 89- level collisional-radiative 
model (CRM). The rate coefficients are obtained by averaging the collision strengths over 
a Maxwellian velocity distribution prevalent in most astrophysical plasmas, unlike the 
approximately monoenergetic velocity distribution in EBITs. Although our line ratios 
are in good agreement with EBIT experiments, other atomic processes, particularly 
cascades from higher levels, may be significant in astrophysical environments. 

Fig. 14 is an example of the 89-level CRM results with the present relativistic 
excitation rates for Fe XVII , and their potential applications to X-ray astronomy in 
sources such as stellar coronae. The detailed collision strengths for the X-ray line 
3F (1-5), an intercombination transition, is shown in Fig. 10(a), with huge resonance 
effects. The temperature dependence of 3F/3C ratio (the X-ray line 3C is an allowed 
transition) is demonstrated in Fig. 14, and compared with previous calculations (filled 
squares). The electron density dependence is small; solid-line and dot-line correspond 
to 10 13 and 10 9 cm respectively. The 4 open circles with error bars are observed 
and experimental values. At all temperatures T < 10 7 K the present line ratio departs 
considerably from previous calculations, to more than a factor of 3 at about 10 6 K — a 
fact of considerable importance in photoionized X-ray plasmas that have temperatures 
of maximum abundance much lower than that in coronal equilibrium T m ~4x 10 6 K 
for Fe XVII , as marked. 

We also find that using our Maxwellian-averaged rates the line ratios are lower 
than some of the astronomical observation data, in particular in the non-flaring active 
region of solar corona fl25| , |5|]. Our theoretical calculations thus support the conclusion 
that some other physical processes contribute to level population kinematics, e.g. level- 
specific recombination-cascades from Fe XVIII to Fe XVII levels, blending of satellite 
lines due to inner-shell excitation of Fe XVI, and resonant scattering of the 3C line 



TT| , |12]| . Investigation of these mechanisms is still underway, however the present 
accurate ab initio theoretical calculations set denifite limits on atomic/astrophysical 
models to estimate opacity effects, column densities, emission measures, and other 
astrophysical quantities. 

7. Conclusion 



The principal features of the present work follows. 

(I) The hitherto most detailed sets of 37CC and 89CC BPRM calculations show that 
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the n— 4 complex explicitly included in the latter calculation has a considerable effect 
on the collision strengths for Fe XVII transitions. In particular, prominent resonances 
appear in the 89CC calculations at energies above and below the n— 3 thresholds, and 
the effective collision strengths are considerably enhanced relative to the smaller 37CC 
calculation. The background collision strengths for some transitions are also affected 
due to inter-channel coupling and consequent re-distribution of flux among the larger 
number of channels in the 89CC case. 

(II) New calculations of atomic structure and transition probabilities for all 89 levels 
have also been carried out, including E1,E2,E3 and Ml and M2 multipole transitions. 
Some results are presented, although the primary focus is on electron excitation. Owing 
to the complexity of Fe XVII relativistic and correlation effects are equally important. 
We note that neither pure LS-coupling nor pure jj-coupling is appropriate for some 
transitions. It is found that the M2 and E3 transition probabilities are sufficiently large 
and may have non-negligible effect on the intensities of some important Fe XVII lines. 

(III) The calculated effective collision strengths are benchmarked against EBIT 
experimental data, and show very good agreement to ~10%, or within experimental 
uncertainties. Resonance enhancements in the intercombination lines 3D and 3E is 
much larger than for the dipole-allowed line 3C, and is crucial to spectral formation 
of these important lines. The strong dependence of line ratios 3C/3D and 3E/3C on 
electron beam energy is explained by the present results. 
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Figure 1. Grotrian energy diagram and important Fe XVII X/UV lines, (a) 
Schematic X-ray transitions to the ground level and X/UV-ray transitions between 
excited levels; (b) Wavelengths and level designations for some X-rays of considerable 
importance in astrophysics and for important soft X-ray laser lines in laboratory 
plasmas are sketched. 
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Figure 2. Collision strength £1 from 37CC BPRM calculations with detailed resonance 
structures as a function of incident electron energy Ef. (a) magnetic quadrupole M2 
transition 1-2; (b) monopole transition 1-15. The green dots and sqaures are RDW 
values in |l8j and in [|l9|, respectively. 
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Figure 3. Partial wave (PW) collision strength fij and partial sum of PW collision 
strength Yli/2 f° r El transitions, (a) and (b): transition 3C(l-27); (c) and (d): 
3G(l-3). For each transition, £lj and are calculated at four Ef. 100 Ry (black 

curve), 200 Ry (green), 300 Ry (blue), and 400 Ry (Red). The filled circles in (b) and 
(d) are the corresponding total collision strengths that partial sum of PW collision 
strength should converge onto. 
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Figure 4. Partial wave (PW) collision strength fij and partial sum of PW collision 
strength ^ or -^1 transition between excited levels 2-44. (a) and (b): from 

89CC; (c) and (d): Vlg and J2i from RDW Q. Other symbols are the same as 
Figs. 3. 
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Figure 5. Partial wave (PW) collision strength flj and partial sum of PW collision 
strength Yli/2^J f° r ^2 transitions between excited levels 2-51 and 6-60. Other 
symbols are the same as Figs. 3. 
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Figure 6. Reduced collision strengths f2 r are plotted as a function of reduced electron 
energy E r = Ei/(Ei + c) (c=e is a constant, E r £ [0,1]) for E2 transitions 2-51 and 
6-60. The infinite energy red filled circles and our RDW blue squares are also shown 
for comparisons. 
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Figure 7. Comparisons of collision strength £1 with detailed resonance structures 
versus incident electron energy Ei between 37CC and 89CC BPRM calculations for 
dipole allowed El transition 3C (1-27, spin unchange). (a) 37CC; (b) 89CC; (c) 
37CC (red curve) vs. 89CC. In (a) and (b), the green and red dashcd-lines are the 
numerical averaged (NA) and gaussian averaged (GA) collision strengths, respectively; 
the green dots, filled squares, open squares and filled and open blue triangles are RDW 
values iQ , RDW values |l9) , DW values [Q and previous R-matrix values g and , 
respectively. The arrows in (c) represent 37CC and 89CC thresholds. 
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Figure 8. Comparisons of collision strength with detailed resonance structures 
versus incident electron energy Ei between 37CC and 89CC BPRM calculations for 
intercombination transition (dipole allowed spin-changed) El transition 3D (1-23). The 
other symbols are the same as Figs. 7. 
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Figure 9. Comparisons of collision strength with detailed resonance structures 
versus incident electron energy Ei between 37CC and 89CC BPRM calculations for 
intercombination transition (dipolc allowed spin-changed) El transition 3E (1-17). The 
other symbols are the same as Figs. 7. 




Figure 10. Collision strength fl from 89CC BPRM calculations with detailed 
resonance structures as a function of incident electron energy Ei for pure El transitions 
from the ground state: (a) 1-5; (b) 1-33. The other symbols are the same as Figs. 7. 



33 




CO 



CO 



CM 




CO 



to 



no 




CO 



to 




co cc 



to c: 



cm 



in 



CO 



to 




CO = 




CO 



CO 



CO 



to 



CM 



CO 



CM 



^5 m6u©Jis uojSjHOQ pe^rei&Q 



Figure 11. Collision strength from 89CC BPRM calculations with detailed 
resonance structures as a function of incident electron energy Ei for E2 and Ml 
transitions: (a) pure E2 (1-7); (b) pure E2 (1-37); (c) pure Ml (1-6); (d) pure Ml 
(1-9); (e) mixed E2+M1 (2-3); (f) mixed E2+M1 (2-5). The other symbols arc the 
same as Figs. 7. 
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Figure 12. Collision strength from 89CC BPRM calculations with detailed 
resonance structures as a function of incident electron energy Ei for El, M2 and E3 
transitions: (a) pure M2 (1-2); (b) pure E3 (1-20); (c) mixed M2+E3 (3-36); (d) mixed 
E1+M2+E3 (2-35). The other symbols are the same as Figs. 7. 
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Figure 13. Collision strength fl from 89CC BPRM calculations with detailed 
resonance structures as a function of incident electron energy Ei for monopolc 
transitions from the ground state: (a) 1-15; (b) 1-11. The other symbols are the 
same as Figs. 7. 
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Figure 14. Intensity ratios of two X-ray lines 3F/3C (1-5/1-27) from Fe XVII as a 
function of electron temperature, calculated at two electron densities 10 13 cm~ 3 (solid 
lines) and 10 9 cm~ 3 (dotted lines), and compared with observed and experimental 
values: from the solar corona at T m ~ 4 million degrees Kelvin ]25|], from the corona 
of solar-type binary star Capella at ~ 5 million degrees Kelvin |26| , and from the 
EBIT experiment at ~ 10 million degrees Kelvin |p^| . The filled squares are values 
using previous cross sections jjj which differ from observations at low temperatures. 



